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Abstract 

The analytical ansatz was found for Bogolubov's angle 9(p) in the 
't Hooft model at N — > oo. The appropriate calculations were done and 
final numeric approximation was found for this angle 9{p). 
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1. Introduction 



There is a well known generalization of QCD from three colors to N (N G N) colors 
or, in other words, from an SU (3) gauge group to an SU (N) gauge group which 
was done by t'Hooft t' Hooft made this generalization because he hoped that 
model with a large N (strictly speaking with limit N —>■ oo, the so-called t'Hooft 
limit) might have an exact solution and might be qualitatively and quantitatively 
close to the actual QCD with N = 3. His hope was justified; theory in the 
N —>■ oo limit has a considerable simplification. The fact is that we should take 
into account only planar diagrams, because others vanish [2] : "In the actual world 
quarks belong to the fundamental representation of SU(3). If we assume that this 
assignment stays intact in multicolor QCD, each extra quark loop is suppressed 
by 1/N. Therefore, in the 't Hooft limit each process is saturated by contributions 
with the minimal possible number of quark loops". This is the main reason we are 
motivated to investigate this model and to find a solution for Bogolubov's angle. 
All definitions and formulas which were necessary for presented calculations were 
taken from the paper [2]. 



2. Analytical ansatz for 9(p) 

The starter formulas for further calculations: 

E p cos6(p) — m = — J dkcos9(k)- 



2 J (p- k) 
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E p sin 6(p) — p = — J dk sin 6(k)- 



2 J v J (p-ky 

with the boundary conditions 
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— at p — > oo 

0{P)^{ I (2) 
— — at p — » — oo 

determined by the free-quark limit. This set of equations was firstly obtained by 
Bars and Green [Hj. The angle 9{p) is referred to as the Bogolubov's angle, or 
more commonly, the chiral angle. 

From these equations one can easily get the following integral equation for the 
Bogolubov's angle OEj: 

7 f 1 

p cos 9(p) — m sin 9(p) = — / dk sin {9{p) — 9{k)) — - . (3) 

2 J [p — k) 

Assuming that the chiral angle was found the following integral equation for the 
E p was obtained in paper [3J E] : 

if 1 

E p = m cos 9(p) +psmO(p) H — / dk cos (9(p) — 9(k)) — — . (4) 

2 J (p-ky 



2 



All calculations were performed in the limit of quark mass m — 0. Thus, one can 
rewrite Eq. (jHJ) and Eq. (0J) in this limit. 



E p = psin6(p) + — 



7 f cos (9{p) — 9{k)) 



dk . 



P 1 2 J {p-kf 

The exact singular analytical solution of the integral equation © is 
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where sign(p) is the sign function 

sign(p) = #(p) — p) . 



(5) 
(6) 

(7) 
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Figure 1: The analytical solution for Bogolubov's angle 8(p) vs. p. 
Substituting this solution in Eq. (JHJ) we can get 



E n 



However, this analytical solution is unphysical [2] for several reasons. For example, 
E{p) becomes negative at \p\ < ^Ff. This feature of the solution (JHJ) cannot be 
amended by a change in the infrared regularization. In fact, solution (jHJ) does not 
correspond to the minimum of the vacuum energy [SJ. 

A stable solution has the form depicted in Fig. 2. It is smooth everywhere. At 
\p\ <C y/j it is linear in p. 

For the convenience of the further calculations, the system of units in which 7 = 1 
was used. The asymptotic behavior of the physical 8(p) [2]: 

6{p) ~ p for \p\ < 1 aM 

-3= (3) 3 + --- for \p\ >1 . ' * <)! 



9{p) 



I sign (p) - -Z= 
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0(p) 
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Figure 2: A stable solution for Bogolubov's angle 9{p) vs. p. 



The asymptotic behavior of 9(p) at \p\ ^> 1 was calculated by calculating the 
chiral quark condensate (^pip) for the smooth physical solution as a self-consistency 
condition [S]. 

Taking these asymptotics into account, we can find the analytical function for 9(p) 
which we can use to fit physical 6(p). Therefore the analytical ansatz for 9(p) is 



m \ ® a P 

9{p) = 7 =-= — - + arctan (ap) 

2 + vrv^a 3 + 6a 2 p 2 



(10) 



where a is a free parameter and a > 0. 

For any a, this ansatz gives the exact asymptotic behavior of the physical 9{p): 



9{p) 



1 + 



6 



2 + vrv^a 3 



ap + o 1 + 



108 



(2 + 7iV6a 3 )' 



a 3p3 _| for |p| <C 1 



9{p) 



and (11) 



7T 1 Svra 3 (3vra 3 + 2^) - 8 i 



The graph of the physical 9(p) from this ansatz looks exactly like the curve 
on Fie. 2. 



3. Numeric solution of 0(p) 

Mathematica 5.1 was employed to create a computer program for numerical calcu- 
lations. Also for these calculations the analytical ansatz was used and it was found 
that, within precision of the calculations, the best fit of the analytical ansatz for 
9(p) is reached when a = 0.3701. 
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The precision of the calculations is 0.0001. The discrepancy of a solution of the 
integral equation (JHJ) at 7 = 1 is calculated by the formula: 



S(a,p) = 

- 1 t /y (r>. COH , f„. a ) - 1 ? sin^ calc fa,q) - 6> calc (fc,a)1 ,, V (12) 

where n is a number of points in which 9 ca \ c (pi) was calculated. 
The dependence of S(a, p) on momentum p is an indirect dependence. This depen- 
dence is realized by the dependence of S(a,p) on the choice of range of momentum 
p and on the fragmentation of this range. 
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Figure 3: The dependence of the solution discrepancy of the integral equation for 
6(p) on the parameter a. 



The calculations of the solution discrepancy S(a,p) shows that there are two lo- 
cal maximums of this discrepancy dependent on momentum p near p = 0.4 and 
p = 1.7 momentum units. Investigating the dependence of the discrepancy on pa- 
rameter a at this region of momentum p, it was found that the discrepancy S(a,p) 
of the calculations is bounded above 0.035, and that it reaches the minimum at 
a = 0.3701. 

For the case a = 0.3701, the solution for the physical Bogolubov's angle 6{p) was 
obtained: 



9(p) = p + arctan (0.3701») 

KFJ 2.3901 + 0.8218p 2i ^ v FJ 



(13) 
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Figure 4: The best fit at a = 0.3701 for the physical solution of the Bogolubov's 
angle 9{p) vs. p. 



The asymptotic behavior of physical 9(p): 
9(p) = 1.2992p - 0.3364p 3 + • • • 
9{p) = I sign(p) 5.9502^ + 

4. Conclusions 

The essential advantage of the suggested technique for calculations of the 9(p) is 
that they are based on the exact asymptotic behavior of the physical Bogolubov's 
angle 9{p) at p — * and at p — > ±oo. This gives us an opportunity to call the 
obtained solution for a semi-analytical one. 

The accuracy of my calculations is 0.035 or better. This error appears when 
momentum p is about either 0.4 or 1.7 momentum units. However, when the 
absolute value of momentum is less than 0.02 or greater than 4.9 momentum 
units, the accuracy has already been an order of magnitude better, i.e. an accuracy 
about 0.003. 

If less of the solution discrepancy of the integral equation for 9{p) is needed, we 
suppose the accuracy of calculations can be improved by adding terms of high 
order of p to the trial function of 9(p). 



for \p\ 1 
for \p\ ^> 1 . 



and 



(14) 
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